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2^ ' Abstract. Our understanding of the cosmic history of galaxy clusters 

^_^ , has recently been enhanced due to an extensive series of observations in- 

^_, • eluding faint spectroscopic data (especially those obtained at the Keck 

O , Observatory), deep optical and NIR imaging from the ground and in 

^J ' space, morphological data from HST, and new constraints on the evo- 

vQ , lution of the intracluster medium from ROSAT and ASCA. When such 

psj I observations are applied to complete, objectively derived catalogs of clus- 

ters, our constraints on cluster formation and evolution become quite con- 
fined. A picture is emerging in which the bulk of cluster formation starts 
at z > 2, yielding cluster potentials that are well established by 2 ~ 1, 
r^ I and there is little substantial evolution of the cluster galaxy population 

Tij- . since z ~ 0.4. This review talk will summarize the current observational 

^D I constraints on the properties and evolution of high redshift clusters and 

11/ ■ protoclusters and their implications. 

5' 

^: 

' , 1. Introduction 

H ; 

'^ I The evolution of galaxy clusters is quite sensitive to the physical processes which 

d ■ dominate the formation of structure and to the cosmological parameters. Cluster 

^ I evolution is inherently complex both because clusters are not closed systems and 

V~j ■ because the 3 main mass components (dark matter, intracluster gas, and galax- 

r> I ies) evolve differently. As a consequence, different cluster parameters evolve on 

C^ ■ different timescales depending on the thermal and dissipative properties of the 

mass component (s) which most strongly control each cluster parameter. Fur- 
thermore, it is becoming evident that the process of cluster formation extends 
over a moderately broad range in redshift. Given the complexity of the prob- 
lem, breakthroughs in understanding the formation and evolution of clusters of 
galaxies must rely on observations spread over a large range of wavelength and 
redshift. 

With the advent of the Keck 10m telescopes, the restored resolution of 
HST imagery, improvements in IR arrays, and the enhanced x-ray imaging and 
spectroscopic capabilities of ROSAT and ASCA, constraints on the properties 
of the ICM and the cluster galaxy population have been extended out to 2: ~ 1 
and beyond. Indeed, cluster candidates have been identified out to z ~ 3 (see 
Table 1). The large temporal baseline these data cover now allow much tighter 
constraints on scenarios for cluster evolution. 



2. Search Strategies 

The most distant cluster candidates have been found by conducting searches in 
the vicinity of high-z radio galaxies or quasars. The evidence for the presence 
of a cluster in these cases is typically based on fewer than 5 spectroscopically 
confirmed members plus a statistical excess of red galaxies or Lya emitters. 

Table 1. Examples of z > 1 Cluster Candidates 



Name 


Redshift 


Reference 


MRC 0316-257 


3.14 


Le Fevre et al. 1996 




QSO 0953+545 


2.50 


Malkan et al. 1996 




QSO 1312+4237 


2.50 


Campos et al. 1999 




53W002 


2.39 


Pascarelle et al. 1995 




QSO 2139-4434 


2.38 


Francis et al. 1996 




3C294 


1.79 


Dickinson et al. 1999 




RXJ0848+4453'' 


1.27 


Stanford et al. 1997; Rosati et al. 


1999 


3C324 


1.21 


Dickinson 1997 




AXJ2019+112 


1.01 


Benitez et al. 1998 




3C184" 


1.00 


Deltorn et al. 1997 





" Spectroscopic confirmation based on more than 10 z's 

Hall & Green (1998) have also performed a search around a sample of radio- 
loud quasars and have identified 31 possible clusters in the range 1 < z < 2. The 
physical properties of z > 1 cluster candidates are not well quantified because 
of the limited amount of spectroscopic data available to date. While looking 
for clusters in the vicinities of radio galaxies or quasars is fruitful, the resulting 
samples will naturally suffer from selection effects associated with limiting ones 
search to such interesting environments. None the less, it appears that over- 
densities which may be the progenitors of present day clusters exist at z ~ 3 
(~ 15% of the current age of the universe). 

A more complete picture of the properties of high-redshift clusters can be 
obtained at z < 1.3 through objective searches of wide areas of sky in optical, 
NIR, and x-ray passbands. The advantage of x-ray cluster selection is two-fold: 
1) emission from the hot intracluster medium (ICM) directly indicates the pres- 
ence of a gravitationally bound system and 2) the ICM comprises 70 to 80% 
of the cluster's baryonic mass. Nearly all x-ray selected high-z clusters are rich 
and elliptical dominated. The Extended Medium Sensitivity Survey (EMSS; 
Gioia et al. 1990a; Henry et al. 1992) has been used to identify clusters out to 
z ~ 0.85 over an area of ~ 850 deg^ and the ROSAT Distant Cluster Survey 
(RDCS; Rosati et al. 1998) has been used to find systems out to z ^ 1.3 over 
an area of ~ 30 deg^. The RDCS, in particular, includes 100 spectroscopically 
confirmed clusters. Of these, 33% have z > 0.4 and 25% have z > 0.5 (Rosati 
1998). While past x-ray telescopes have had fairly low effective areas, new ob- 
servatories, like XMM, will provide at least an order of magnitude improvement. 
Equally exciting are developments in the use of the Sunyaev-Zeldovich Effect to 
locate clusters. Mohr et al. (1999) indicate that SZE facilities in the near future 
will be able to detect ~ 100 z > 1 clusters per year. 

Searching for distant clusters in the optical and NIR, however, also has 
significant advantages. From a practical point of view, there are more telescopes 



and larger area mosaic cameras available in the optical/NIR than in x-rays. 
Optical/NIR searches will also find clusters spanning a wider range of x-ray 
luminosity and total mass {e.g., Holden et al. 1997). Although the spurious 
detection rate at high-z can be ~ 30%, the use of photometric redshifts can 
dramatically reduce the number of false positives. Some of the largest area 
and deepest optical/NIR distant cluster surveys include the Palomar Distant 
Cluster Survey (5.1 deg^; Postman et al. 1996), the ESO Imaging Survey (12 
deg^; Scodeggio et al. 1999) The Deeprange Survey (16 deg^; Postman et al. 
1998), and the NOAO Deep- Wide Survey (18 deg^; Jannuzi & Dey 1999). 

An optimal strategy, of course, is to combine x-ray and optical/NIR data 
obtained over the same region of sky. This allows a full assessment of the se- 
lection biases to be made and is likely to reveal subtle effects which can be 
important in interpreting, for example, the abundance of high-z clusters. The 
benefits of such joint searches are already being realized: Scharf et al. (1999) 
used 22 deep ROSAT PSPC fields as targets for deep optical imaging to study 
the effects of optical and x-ray selection on derived cluster evolution and to look 
for correlations in the large-scale distribution of diffuse x-ray emission and the 
galaxy distribution. Preliminary results include the possible first direct detec- 
tion of x-ray emission from an intercluster filament at z ^-^ 0.4 — 0.5. Stanford et 
al (1997) and Rosati et al (1999) have identified a supercluster at z = 1.27 in the 
Lynx field which was initially detected in the NIR (K-band) and, subsequently, 
in x-rays. 



3. Cluster Abundance 

The abundance of clusters as a function of redshift is one of the fundamental 
constraints on both structure formation and cosmological models. The space 
density of clusters at z > 0.5, for example, is highly sensitive to 0,^ (Bahcall, 
Fan, & Renyue 1997; Donahue & Voit 1999). Present observational constraints 
from x-ray surveys (followed up by optical spectroscopy) indicate that the co- 
moving space density of clusters per unit Lx is invariant out to at least z = 0.8 
for systems with L^ < 3 x 10^^ erg sec^^ (Henry et al. 1992; Rosati et al. 
1998). For more luminous (massive) clusters, mild negative evolution has been 
reported (Gioia et al. 1990b; Henry et al. 1992; Vikhlinin et al. 1998) although 
the deficit, expressed in absolute numbers, is only a dozen or so EMSS clusters 
at z > 0.4 (small enough that one might worry about subtle selection biases at 
low x-ray surface brightness levels in the existing surveys). The distribution of 
poor to moderately rich optically selected clusters is also consistent with a con- 
stant comoving space density to at least z = 0.6 (Postman et al. 1996; Holden 
et al. 1998). 

At z > 1, our constraints on cluster abundances presently suffer from a 
lack of data. There are at least 5 known clusters with 0.75 < z < 1.3 that 
have velocity dispersions in the range 700 < aiD < 1400 km sec^^. At least 
two of these, MS-1137 {z = 0.78) and MS-1054 {z = 0.83), have relatively high 
kinetic gas temperatures (Tx) - 5.7 keV and 12.4 keV, respectively (Donahue et 
al. 1998, 1999). The existence of massive (> 5 x 10^"^ M©) clusters at z ~ 1 
is, thus, no longer in doubt. Interestingly, the space density of z ~ 1 clusters 
inferred from the RDCS coverage of the Lynx region is '-^ 2.4 x 10~^/i^ Mpc~^ 
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Figure 1. The Donahue et al. (1998) T^ - oxn relationship for 6 rich 
clusters in the range 0.5 < z < 0.83. The best fit low-z relationships 
are shown for comparison. The distant cluster data are consistent with 
that seen at z < 0.2. 



(Rosati 1999), within a factor of 2 of the density z ~ 3 structures delineated 
by Lyman break galaxies (Steidel et al. 1998). Whether this is indicative of an 
evolutionary connection or mere coincidence remains to be decided. 



4. Evolution of the Gravitational Potentials of Clusters 

The evolution of the ICM and its correlation with global cluster kinematics 
provide direct constraints on the growth of the gravitational potentials which we 
call clusters. The relationship between L^ and T^ in low-z clusters is remarkably 
tight but is somewhat steeper than that predicted by bremsstrahlung emission 
from a population of virialized, structurally identical clusters with constant gas 
fraction (Arnaud &: Evrard 1999). However, with reasonable constraints on 
cluster structure, the same authors find the fractional variation in cluster gas 
fraction is < 15%. Mushotsky & Scharf (1997) demonstrate that the L^ — T^ 
relationship exhibits no significant evolution out to z ~ 0.4. Donahue et al. 



(1999) have extended this work out to z ~ 0.9: if the evolution of the relation 
is parameterized as L^ = T^{1 + z) , then they find that A > 1.5 is rejected 
with greater than 3a confidence (for Qq = 0.5). Values of A = 2 — 3 would be 
required to explain the lack of evolution in the x-ray luminosity function cited 
above if Qm = 1- Donahue et al. (1998) have further shown (see Figure 1^) that 
the relation between the cluster velocity dispersion {(Tid) and T^ is invariant 
out to z ~ 0.8. Cluster potentials are clearly well established in the universe by 
z ~ 0.9 and, on average, the x-ray properties of the ICM are similar to those in 
current epoch clusters. 

The distribution of the optically luminous mass in clusters, as delineated 
by the member galaxies, may be experiencing more recent evolution than the 
ICM. Clusters exhibit significantly more asymmetry in their galaxy distribution 
at z > 0.7 than at the present epoch (Lubin & Postman 1996) - the observed 
profiles are inconsistent with azimuthal symmetry at the 99.9% confidence level, 
in strong contrast with the situation at z < 0.3. In some cases, like MS-1054, 
the dumpiness seen in the galaxy distribution is also seen in the x-ray brightness 
distribution (Donahue et al. 1998) and in a mass map based upon weak lensing 
distortions (Hoekstra, Franx, & Kuijken 1999) - characteristic of recent merger 
activity. Indeed, mergers of group-size clumps at z ~ 1 may be the origin of 
some of the current epoch richness class 1 and clusters (Lubin et al. 1998; 
Gioia et al. 1999). However, the majority of the known high-z clusters appear 
to have been in existence since at least z ~ 2, as observations discussed below 
suggest. 



5. Constraining the Epoch of Cluster Formation 

The processes which control the formation of clusters leave observable signa- 
tures in the evolution of the morphological and spectrophotometric properties 
of cluster galaxies. This is a key reason why observational work in this area 
has been a major component of recent extragalactic telescope programs. The 
evolution of the mass function of cluster galaxies, in particular, provides crit- 
ical constraints on and tests of cluster formation scenarios. At high-z, cluster 
galaxy mass determinations are difficult to obtain but the K-band cluster galaxy 
luminosity function (KLF) can provide a reliable substitute because it probes 
the total stellar mass component and is not strongly sensitive to the instanta- 
neous star formation rate (e.g., see Gavazzi, Pierini, Boselli 1996). De Propris 
et al. (1999) have derived the KLF for 38 clusters at 0.1 < z < 1. Their main 
result is that the KLF departs from no-evolution predictions at z > 0.4 how- 
ever the changes observed are consistent with simple, passive evolution (aging 
of the existing stellar population) and a narrow formation epoch around z = 2 
(if A = 0.7) or z = 3 (if A = 0). 

Comparison of the broadband colors and spectral features of the early type 
cluster members at 0.76 < z < 0.92 with spectral synthesis models suggests these 
galaxies are old (mean ages ~ 3 it 2 Gyr, at the observed redshifts) implying 
a relatively early formation at z > 2 as well (e.g.. Bower, Kodama, Terlevich 
1998; Postman, Lubin, Oke 1998; Stanford, Eisenhardt, Dickinson 1998). Such 
observations also suggest that the mass-to-light ratios of the early type cluster 
galaxies have evolved passively since at least z ~ 1.2 (Kelson et al. 1997; van 
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Dokkum et al. 1998). Taken in concert with the results on the KLF evolution, 
one may conclude that the mass function of cluster galaxies has remained roughly 
invariant since z ~ 1.2. 

An additional constraint of the duration of the cluster galaxy formation 
era comes from the optical/NIR color-magnitude relations for the red galaxy 
population which are well-defined and exhibit remarkably low scatter in clusters 
from < z < 1 (Stanford, Eisenhardt, Dickinson 1998). This places a stringent 
constraint on their formation synchronicity of Ai < 4 Gyr (roughly the time 
between 2; = 10 and z = 1.5 in a Oo = 0.2, h = 0.6, A = cosmology). The 
coeval nature of cluster elliptical evolution is also reflected in observations which 
are consistent with exponentially decaying star formation rates with relatively 
short e-folding times (0.1 < r < 0.6 Gyr; Postman, Lubin, Oke 1998). 

There is some evidence (Fuller, West, Bridges 1999) that the brightest clus- 
ter galaxies are preferentially aligned with the global cluster galaxy distribution, 
an effect also suggestive of an early formation epoch. 

The metallicities of cluster ellipticals in the redshift range 0.5 < z < 1 are 
consistent, on average, with the close to solar values observed in current epoch 
ellipticals. Similarly, there does not appear to have been much change in the 
metallicity of the ICM (~ 0.2 — 0.45 solar) between now and z ~ 0.8 (Donahue 
et al. 1999). This suggests reprocessing of the baryonic mass component of 
clusters had been on-going for a few Gyr prior to current lookback time. Hier- 
archical models do predict, however, that some star formation activity should 
be occurring in clusters since z ~ 1. This is indeed seen in the spiral members 
and in the surrounding field galaxy population. For example, the fraction of 
cluster members with strong Oil emission (EW > 15A), a reliable star forma- 
tion indicator, increases by a factor of 3 - 4 between now and z ~ 0.92 (Balogh 
et al. 1998; Postman, Lubin, Oke 1998). The percentage of cluster galaxies 
with post-starburst spectral features increases nearly tenfold between now and 
0.3< z < 0.5 (Dressier et al. 1999). 



6. Evolution of the Morphological Mix 

Using the high angular resolution imaging provided by the Hubble Space Tele- 
scope, several teams (e.g.. Dressier et al. 1997; Oemler et al. 1997; Small et 
al. 1997; Lubin et al. 1999) have conducted morphological surveys of clusters 
in from 0.3 < z < 1. Upon comparison with similar ground-based studies of 
low-z clusters, it appears that the distribution of cluster galaxy morphologies 
has undergone rather substantial evolution between 0.4 < z < 1 but remains 
relatively invariant between z = 0.4 and the present epoch. One indicator that 
has been used to gauge the evolution is the ratio of the number of lenticulars 
(SO) to ellipticals in the central regions of clusters. Figure ^ summarizes the 
current constraints on the redshift dependence of the SO/E ratio in clusters ob- 
served by Dressier et al. 1997 (D97) and Lubin et al. 1999. Although the D97 
results have been widely cited as evidence for substantial morphological evolu- 
tion even since z = 0.4, re-analysis of the data by Andreon (1998) and new data 
at higher z from Lubin et al. 1999 suggest that the metamorphosis may not be 
as dramatic as originally thought. However, the high-z data from Lubin et al. 
also find that the while ellipticals still preferentially reside in the highest density 
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Figure 2. The evolution of the SO/E ratio based on WFPC and 
WFPC2 imaging of clusters in the range 0.3 < z < 0.95. The low- 
z reference is from Dressier (1980). The intermediate redshift results 
are from Dressier et al. 1997. Andreon (1998) re-analyzed their data 
and finds higher SO/E ratios which are consistent with the local value 
- his ratio for CL0939-I-47 is highlighted. The three high-z clusters are 
at z = 0.76, 0.90, and 0.92 (Lubin et al. 1999). 



environments, the Spiral/SO morphology-density relation at z > 0.5 is much less 
well-defined than it is now. 

There is further evidence that morphological modifications are occurring 
at least as recently as z ~ 0.3 — 0.4: Dressier et al. (1999) find that the Oil 
equivalent widths, for a given morphological type, are lower in 0.3 < z < 0.5 
cluster galaxies than in the field and that the actively star-forming galaxies in 
these clusters have a more extended spatial distribution than the non-active 
galaxies. The kinematic properties of the early and late type cluster galaxies 
appear to differ as well, although these differences exist both in the present 
epoch and at z ~ 1 (Adami et al. 1998; Lubin et al. 1999). Specifically, the 
spiral population tends to have a higher velocity dispersion than the elliptical 
members suggestive of an on-going spiral infall process. Poggianti et al. (1999) 
find, however, evidence for widespread cessation of star formation activity in 
intermediate z clusters over a relatively short (~ 1 Gyr) timescale. Specifically, 
90% of the spiral cluster members they studied show spectral signatures of either 
enhanced or suppressed star formation relative to local spirals. 



7. A Possible Scenario 

I will now propose a possible evolutionary sequence which incorporates the myr- 
iad of constraints derived from observations of clusters. Note that some of the 
steps below have not been observationally confirmed or are still controversial! 

• In a universe with 0^ < 1 (and possibly A 7^ 0), protoclusters form at 
z>3. 

• The sites of formation are located at the intersections of filamentary mat- 
ter fiows and the first cluster galaxies form during the first generation of 
matter crossings {This is a conjecture based solely on N-hody simulations 
and popular dark matter models). 

• The richest, current epoch clusters formed first. Some of the poorer clus- 
ters seen today may have developed via group-group mergers since z ~ 1. 

• Primordial ICM shocks and begins to emit x-rays at z ~ 2 (and perhaps 
earlier). Enrichment of the ICM most likely occurs in the z > 1 era. From 
z ~ 1 to now, there is little evolution of the ICM. 

• The brightest cluster galaxies grow via cannibalism until z ~ 1.5. Most 
merger activity ceases by z ~ 1 and subsequent evolution is passive. Other 
massive ellipticals assemble prior to z ~ 2 and are the first to reach dy- 
namical equilibrium with the cluster potential. 

• The most active periods of star formation within the cluster occur at z > 1. 
Most star formation is quenched, however, by z ~ 0.4 — 0.5. 

• Infall of spirals results in morphological and color gradients within the 
cluster. This process continues up to the present epoch. 
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• The SO and dwarf elliptical populations develop within the cluster core, 
certainly by z ~ 0.5 and more likely by z ~ 1. The likely relevant pro- 
cesses involved are ram pressure stripping, mergers, and tidal stress. SO's 
may descend from high surface brightness spirals, dE's from low surface 
brightness spirals (Moore at al. 1998). 

There are notable exceptions to the above scenario such as low-z spiral rich 
clusters (e.g., Virgo, Hercules) and low-z irregular clusters (e.g., Abell 1185), 
which are probably still dynamically young, suggesting that some cluster evo- 
lution is still occurring at the present epoch. Furthermore, our knowledge of 
cluster evolution at z > 1 is still quite rudimentary. Thus, while great strides 
have been made, there remain many steps to go before our understanding of 
the cluster formation process is complete. Some of the observational programs 
which will take us farther towards this goal are now, or soon will be, underway. 
These include more complete and larger z > \ cluster samples, more objec- 
tive and precise studies of the z = cluster population (e.g., the Sloan Digital 
Sky Survey, the 2dF survey, the REFLEX survey), improved x-ray observations 
from XMM and Chandra and optical/IR observations with HST (ACS, WF3) 
and SIRTF, extended spectroscopic studies of high-z clusters using the growing 
suite of 8 - 10m ground-based telescopes, construction of mass-selected cluster 
catalogs from SZE surveys, and ultra deep 21 cm searches for protoclusters at 
z > 2 using the Giant Meter-wave Radio Telescope. 

Acknowledgments. I thank the SOC for their generous travel support 
which made attendance of this meeting possible. I also wish to thank Megan 
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of clusters in advance of publication. 
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